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ABSTRACT 

Recent analyses of the stellar stream of the Sagittarius dwarf galaxy have claimed that the kinematics 
and three-dimensional location of the M-giant stars in this structure constrain the dark matter halo of 
our Galaxy to possess a triaxial shape that is extremely flattened, being essentially an oblate ellipsoid 
oriented perpendicular to the Galactic disk. Using a new stream-fitting algorithm, based on a Markov 
Chain Monte Carlo procedure, we investigate whether this claim remains valid if we allow the density 
profile of the Milky Way halo greater freedom. We find stream solutions that fit the leading and 
trailing arms of this structure even in a spherical halo, although this would need a rising Galactic 
rotation curve at large Galactocentric radius. However, the required rotation curve is not ruled out 
by current constraints. It appears therefore that for the Milky Way, halo triaxiality, despite its strong 
theoretical motivation, is not required to explain the Sagittarius stream. This degeneracy between 
triaxiality and the halo density profile suggests that in future endeavors to model this structure, it 
will be advantageous to relax the strict analytic density profiles that have been used to date. 
Subject headings: Galaxy: halo — dark matter — Galaxies: dwarf 



1. INTRODUCTION 

In standard A Cold Dark Matter (A-CDM) cosmology 
(jKomatsu et al.lfeOllf ) galaxies form within dark matter 
halos that coalesce through repeated mergers. High reso- 
lution pure CDM simulations give rise to significantly tri- 
axial galaxy halos, that become progressively less triaxial 
towards lower sizes and masses due to the increased dy- 
nami cal age of the lower mass structures (|Allgood et al.l 
I2QQ6I ). The inclusion of baryons in galaxy formation 
simulations alters dramatically the behavior of the cen- 
tral regions of the halo, which added to the effect of 
merging satellites make s the central halo become rounder 
(|Debattista et al.ll2008f ). 

In the Milky Way, we have perhaps the best oppor- 
tunity to constrain the shape of any dark matter distri- 
bution. One particularly promising strategy to accom- 
plish this, and which is generally not possible in extra- 
galactic systems, is to use streams of stars to probe the 
global potential. The streams of low-mass satellites fol- 
low fairly closely the centre of mass orbit of their pro- 
genitor, and with careful modeling it is possible to find a 
solution for the stream structure within a Galactic mass 
distribution. This approach has been used in the past to 
model the stream of the Sagittarius (Sgr) dwarf galaxy 
(jlbata et al.ll2001f h probably the most significant accre- 
tion into the Milky Way in the last ~ 5 Gyr. Recently, 
lLaw fe M ajewski ()2Q1QL hereafter LM10) presented a 
new analysis of the spatial and kin ematic structure of M- 
giant star members of this stream (M ajewski et al.|[2QQ3l 
120041 ). and showed that their data could be reproduced 
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if the dark halo of the Milky Way has a triaxial form 
with a minor to major axis ratio of (c/a)$ = 0.72 and 
intermediate to major axis ratio of (b/a)<$> = 0.99 (in the 
potential). The structure is therefore almost an oblate 
ellipsoid, but with its minor axis contained within the 
Galactic plane, which presumably induces non-circular 
orbits at radii where the disk is not dominant. As LM10 
point out, such a halo is not natural in A-CDM. It is 
also worth noting that their model has an extremely high 
degree of oblateness in density. These c onclusions have 
been f urther reinforced by the analysis of lDeg fc Widrowl 
(|2012l . hereafter DW12) who in fitting the stream with or- 
bits in a Galactic mass model, found a similar oblate dark 
matter halo also oriented perpendicular to the Galactic 
plane. In this contribution we aim to determine whether 
it is possible to reproduce the structure and kinemat- 
ics of the Sgr stream, at least as well as LM10, without 
invoking a triaxial Galactic halo. 

2. MODELLING 

One means to model the formation of the Sgr stream 
would be to make a mass model of the Milky Way and 
an N-body model of the Sgr dwarf, and to evolve the 
latter inside the potential provided by the former using 
an N-body integrator. Through a judicious sampling of 
orbital and structural parameters of the two bodies it 
would be possible to assess the triaxiality issue. This is 
essentially the strategy adopted by LM10, but it has the 
disadvantage of being computationally extremely costly, 
which greatly limits the resolution of the parameter space 
survey. An alternative approach, adopted by DW12, is 
to follow orbits in a mass distribution, however due to 
self-gravity this is no t a good approximati on for streams 
of massive satellites (Var ghese et al.ll201l[ ). 

Instead, we have chos en to use a new strea m-fitting 
algorithm presented in I Varghese et al.l (j2011f ). which 
searches through the parameter space of the host galaxy 
mass model and the satellite's orbital properties and 
mass using a Markov Chain Monte Carlo (MCMC) 
scheme. As we have shown in that contribution, the 
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Figure 1. In panel (a), the sky positions in Right Ascension and Declination of the M-giant stars studied by LM10 (red circles) are 
plotted on top of the the best-fitting model of the Sgr stream found in 10 6 iterations of the MCMC algorithm (green dots). The model 
points have been convolved with a Gaussian function to mitigate against sampling noise in the model, resulting in the probability density 
function map (black). Regions close to the Galactic plane have been masked out. Panel (b) shows the sky regions covered by the Sloan 
Digital Sky Survey (SDSS), where the stream is detected: our estimates of the SDSS sky positions of the main leading and trailing arms are 
shown with turquoise points. The model points that overlap these area are shown again in green. The lower panels show the corresponding 
distance profile, with (c) displaying the M-giants and (d) the distance estimates to red-clump stars in the SDSS as measured by Correnti 
et al. 2010. The fitted model can be seen to give a good approximation to all these data. The LM10 model (a random 10% selection of 
their particles) is shown in all panels with small purple dots. 



method is able to provide a very close approximation 
to a stream without having to resort to N-body simula- 
tions. The simple insight behind this approximation is 
that a leading stream is formed from stars that leave the 
satellite close to the inner LI Lagrange point, while the 
trailing stars leave close to the L2 Lagrange point. In this 
way, we avoid having to undertake a full N-body simula- 
tion, and can feasibly probe millions of model parameter 
combinations in an automatic way with a MCMC algo- 
rithm. 

Following iDehnen fc Binnevl (|1998l hereafter DB98), 



we model the Milky Way as a sum of density compo- 
nents, using the multipole expansion technique to solve 
the Poisson equation. For the bulge, disk, thick disk and 
interstellar medium component, we have adopted the pa- 
rameter values of models '2b' and '2c' from DB98: these 
are realistic density models that fit a wide range of ob- 
servational constraints. DB98 fit these baryonic com- 
ponents within a halo component modeled as a double 
power law; the main difference between their models '2b' 
and '2c' is that '2b' was forced to have an inner power- 
law exponent 7^ = —2, while model '2c' had 7^ = 1, 
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similar to that of a univers al (NFW) cosmological halo 
model (|Navarro et al.|[T997l ) (naturally, these constraints 
on the halo affect the fitted parameters of the other com- 
ponents). However, instead of adopting the DB98 dou- 
ble power-law dark matter halo model, we allow for a 
much larger amount of flexibility in the halo radial den- 
sity profile, as we describe next. The halo is set to be 
an axisymmetric model, defined by a density at lOkpc 
together with an inner power law exponent 7 inside of 
lOkpc, a density at 60kpc together with a power law 
exponent f3 beyond that radius, and density values at 8 
logarithmically-spaced radial locations between 10 and 
60 kpc. A spline function is used to interpolate the den- 
sity between these anchor points. Since the dynamics of 
the Sgr stream should not be very sensitive to the mass 
distribution in the inner Galaxy, or the distribution in 
the outer Galaxy beyond its apocenter, we decided to fix 
7=1 and f3 = 3, similar to the inner and outer behavior 
of an NFW halo model. However, the 10 halo density 
values are free parameters in our modeling. Although 
our code also allows one to set (or fit) the flattening in 
density (g m ) of the axisymmetric halo component, for 
this particular experiment we chose to probe the case for 
a spherical halo, and so kept q m = 1 fixed. 

The remaining parameters that need to be defined are 
the mass, 3-dimensional distance and velocity of Sgr as 
well as the position and velocity of the Sun, in order to 
correctly project the stream kinematics into line of sight 
observables. 

For the initial anal ysis presented here, we chose 
vsgr,heiio — 141 kms -1 (Bellazzin i et al.l 1200 8). and fol- 
lowing LM10 we adopted ds gr — 28 kpc, and Ms gr — 
2.5 x 1O 8 M0; these values were kept fixed. The to- 
tal space velocity and direction of Sgr was allowed to 
vary, as were the Galactocentric radius of the Sun Rq 
and the circular velocity of t he Local Standard of Rest 
vlsr- To be consistent with M ajewski et al.l (j2QQ4i ). we 
adopted their chosen Solar peculiar velocity of (u, v, w) = 
(-9,12,7) kms- 1 . 

The MCMC algorithm that drives the simulation 
has b een updated from that used in Varg hese et al.l 
(2011]), and now uses an afflne- invariant ensemble- 
sampler scheme ( with parallel tempe ring), identical to 
that described in llbata et al.l ([20121 ). The algorithm 
starts off from an initial guess for the 14 model param- 
eters (10 halo density values, 2 velocity parameters for 
Sgr, as well as R@ and vlsr), to begin its exploration of 
parameter space. Our initial estimations for the MCMC 
steps were adjusted dynamically in order to attain a tar- 
get acceptance ratio of 25%. 

At each iteration, the algorithm makes a stream model 
that consists of a prediction of the locus of a tidal stream 
in space. We convert this to a probability density func- 
tion by convolving the points by Gaussian functions that 
approximate the combined effect of the uncertainties on 
the observations (position, distance, velocity) and the 
expected intrinsic width of the stream. Finally, a likeli- 
hood for this set of parameter values is calculated from 
the product of the probabilities of the data (M-giant po- 
sitions, distances and velocities; and SDSS field positions 
and distances) and the likelihood of t he mo d elled Galac- 
tic rotation given the compilation by [Some (20121). 

3. BEST-FIT MODEL 



The best-fit model in 10 6 iterations of the above algo- 
rithm adopting the bulge, disk, thick disk and interstel- 
lar medium parameters of model '2c' of DB98 is shown 
in Figures HKl Figure [T] displays the sky positions (top 
panels) and heliocentric distances (bottom panels) of the 
MCMC solution. The locus of the stream in these pa- 
rameters is displayed in green, these are convolved with 
the parameter uncertainties and intrinsic stream width 
estimate to produce the probability density function dis- 
played in the background grayscale image. For compar- 
ison we also display the LM10 triaxial N-body model 
(small purple dots). The left-hand panels show the po- 
sitions of M-giant stars (red circles) fro m LM10 within 
15° o f the orbital plane of the Sgr dwarf ([Maiewski et al.l 
120031 ); RA ranges close to the disk or close to the main 
body of the Sgr dwarf have been masked out. We dis- 
carded the region immediately surrounding Sgr as we 
wanted to avoid being affected by the dynamical proper- 
ties of the progenitor (e.g. rotation). Panel (b) shows the 
sky positio n of the main detection of the Sgr stream i n 
the SDSS (iBelokurov et al.ll2006l : IKoposov et al.l 120121) . 
as measured by us directly from SPSS data. Panel 
d) sh ows the distance measurements of iCorrenti et al.l 
20101 ). In both right-hand panels we have masked out 
ranges in RA not probed by the data. It can be seen 
that the MCMC solution follows closely the position and 
distance constraints. 

Figure [2] compares the radial velocity profile of the 
model to the M-giant data: the kinematic behavior is 
similar to that of the LM10 model, and reproduces well 
the trend in the M-giant radial velocities. The three- 
dimensional structure of the data and models is shown 
in Figure [3j which gives an edge-on view of the Milky 
Way. As before, the red circles are the selection of 
LM10 M-giant s , with SDSS red-clump distances from 
ICorrenti et al.l ([2010) shown with orange circles. The 
green dots mark the locus of the MCMC stream so- 
lution: the correspondence with the M-giant distances 
is not perfect in the leading arm, but neither is that 
of the LM10 N-body simulation (purple dots). Note 
that the M -giant distance uncert ainties are considerable 
(~ 20% — iMaiewski et al.ll2003D . and there may also be 
a systematic difference with respect to the measurements 
based on the red clump or main sequence. 

A comparison of the relative likelihood of our model 
versus that of LM10 is not very useful at this stage, the 
models are very different, and indeed both models fail to 
repr oduce the data in deta il (such as the stream bifurca- 
tion [BelokuroveF^l][J006), but it can be seen that in po- 
sition, distance and velocity the two models are similar, 
and give a similarly good representation of the youngest 
parts of the leading and trailing arms. 

While we use very similar constraints and assumptions 
to LM10, we have fit the Sgr stream without having to 
invoke the need for triaxiality. However, the price to pay 
for this is a peculiar radial mass distribution in the halo. 
Figure [4] shows the required rotation curve and corre- 
sponding density profile of the fitted spherical halo, again 
using the (fixed) non-halo components of model '2c' from 
DB9 8. The red do ts in panel (a) mark the measurements 
from Sofue (2012). Beyond the Solar Circle there is sub- 
stantial uncertainty in the rotation curve, as the large 
scatter reveals, and beyond 20 kpc there are virtually 
no direct constraints. Expectations from A-CDM would 
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Figure 2. The heliocentric velocities of the LM10 M-giants (red 
circles) are compared to the best-fitting model (green points). The 
correspondence can be seen to be as good as that of LM10. 

suggest that the rotation curve should decline slowly be- 
yond 10 kpc (jKlvpin et al.ll2002l ). in stark contrast to our 
model's requirement. Nevertheless, we note that HI ob- 
servations of our n eighboring galaxy, A ndromeda (black 
dots in panel 'a' — iChemin et al.ll2009j ), show a rotation 
curve rise beyond 25 kpc that is similar to our model. 

While it was not our intention to attempt to fit an ac- 
curate mass model of the inner Galaxy, it is worth noting 
that the surface mass density of the disk and the vertical 
force K z lie within observationally-acceptable ranges (see 
DB98). In addition, the Oort constants for our model 
(which are not identical to those of DB98 since we have 
altered the halo) are A = 13.8 and B = —13.5, close to 
the observed A = 14 ± 0.8, B = -12.4 ± 0.6 (see Ta- 
ble 3 of DW12 for comparison to other solutions). The 
detailed structure of the mass distribution in the inner 
Galaxy turns out to be relatively unimportant in fitting 
the Sgr stream: if we use the non-halo components of 
model '2b' from DB98 (instead of those of '2c' as in Fig- 
ures 1-4) the MCMC algorithm also converges on similar 
solutions that require a rising rotation curve. 

4. DISCUSSION AND CONCLUSIONS 

The analysis we have presented here shows that it is 
possible to reproduce the spatial and kinematic struc- 
ture of the Sgr stream at approximately the same level 
of precision as LM10 without the need for a triaxial halo, 
if the Galactic rotation curve rises significantly in the 
range ~ 20 kpc to ~ 60 kpc. Interestingly, the break 
in the density profile occurs at approximat ely where the 
outer stellar halo component proposed bv I Car olio et al.l 
(|2007l ) becomes dominant. This rotation curve rise of 
course places a substantial amount of matter in the outer 
Galaxy, which has the effect of braking the leading arm, 
bringing its velocity into agreement with that of the M- 
giants, particularly in the region around RA ~ 200°. 
In addition to the MCMC parameter survey discussed 
above, we have also attempted an interactive investiga- 
tion, choosing parameter combinations judiciously; how- 
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Figure 3. An x-z cut through the Milky Way, showing the stream 
location (green dots) according to our best-fit model. The path of 
the centre of mass of the Sgr dwarf is shown in red; the fact that this 
does not line up with the stream position is due to the self-gravity 
of the dwarf galaxy, which we take into account in our modeling, 
but is neglected in the DW12 analysis. The LM10 M-giants are 
shown with red circles, and the big orange dots mark the Correnti 
et al. 2010 'A'-field measurements from the SDSS. The small cross 
at the origin mark the Galactic Centre, and the position of the Sun 
is shown with a black circle at (—8,0) kpc. The numbers marked 
on the circular ticks are the corresponding Right Ascension values 
(in deg) in the plane of the Sgr stream. 

ever, we were unable to find another acceptable solution 
without a strong rotation curve rise. 

It is unclear at present what the structure of galac- 
tic mass dark matter halos actually is in nature, but 
in A-CDM, dark matter density profiles decline as r -2 
only when averaged over a very large number of halos. 
The profile of Fig. HJd is therefore not inconsistent with 
A-CDM, and indeed a rotation curve rise of similar mag- 
nitude is observed in M31 over a similar radial range 
(Fig. 2^). We suspect (although we have not modeled 
this) that there is a trade-off between triaxiality and ro- 
tation curve rise, so that a less triaxial model than that 
of LM10 (which is actually an oblate spheroid oriented 
perpendicular to the Galactic disk) can be fitted with a 
less extreme rotation curve than that presented here. 

A consequence of this new solution is that the Milky 
Way would be considerably more massive than estimated 
from models with a falling rotation curve within R < 
60 kp c. Assuming a virial radius of 258 kpc ([Klvpin et al.l 
12002( 1. leads to a virial mass of 2.6-3.1 x 10 12 M for our 
model. It is interesting to note that the LMC and SMC 
would clearly be bound to the Milky Way if this reflects 
the true mass of our Galaxy, contrary to conclusions 
based on HST proper motion measurements that assume 
a sta ndard Galactic mass model (see. e.g., iBesla et al.l 
120071 ). To be consistent with the timing argument limit 
for the mass of the Local Group of 5. 58±p x 10 12 M 
(|van der Marel fc Guhathakur ta 2008), the dark matter 
would have to be primarily concentrated within the two 
giant galaxies. While the above mass for the Milky Way 
is certainly substantially higher than most recent esti- 
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Figure 4. The rotation curve required by our non-parametric spherical halo model is shown with a solid line in panel (a), with its 
decomposition i nto halo (das hed line), bulge (dot-dashed line) and disk (dotted line) components. This rotation curve fits the available 
measurements of Sofue (2012) (red dots), but it requires significantly rising values beyond ~ 30 kpc. Altho ugh such a risin g rotation curve 
is unexpected, it is similar to the observed HI rotation curve of the Andromeda galaxy (black circles — IChemin" et al.l[2009h . The two error 
bars at Rqc — 18 kpc and 60 kpc show the lcr uncertainties derived from the MCMC algorithm (note that a prior impedes solutions with 
Vdrc > 300kms — 1 ). Panel (b) shows the corresponding density profile (solid line). The error bars mark the radial location of the density 
anchor points together with the lcr spread from the MCMC fit. For reference we also show an NFW profile (dots) of scale radius 21.5 kpc 
(Klvpin et al. 2002): while this has an identical asymptotic behavior to the halo model shown, it is substantially different in the radial 
range probed by the Sgr stream. 



mates (e.g., I Wat kins et al.ll2010l ). almost all the mass 
constraints at large distance are based on Jeans equations 
solutions that make the unjustified assumptions of sta- 
tistical independence of the halo tracers and some form 
for the orbital anisotropy. 

In summary, we present this solution as a means to 
avoid the triaxial halo models of LM10 or DW12, and in 
particular their peculiar orientation as oblate structures 
perpendicular to the Galactic plane. We recognize of 
course that the solution presented here is also strange, 
being spherical and possessing a very high mass, but we 
judge that it is not ruled out by observational constraints. 
The significant shortcomings of both the LM10/DW12 
and our models suggest that they do not include sufficient 
of the actual complexity of (probably both) Sgr and the 
Milky Way halo. Future modeling of the more ancient 
parts of the Sgr stream may help discriminate between 
these models, and pave the way to better solutions. 
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